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Abstract 

Baryon Acoustic Oscillations (BAO) are frozen relics left over from the pre-decoupling 
universe. They are the standard rulers of choice for 21st century cosmology, provid- 
ing distance estimates that are, for the first time, firmly rooted in well-understood, 
linear physics. This review synthesises current understanding regarding all aspects 
of BAO cosmology, from the theoretical and statistical to the observational, and 
includes a map of the future landscape of BAO surveys, both spectroscopic and 
photometriclil. 



1.1 Introduction 

Whilst often phrased in terms of the quest to uncover the nature of dark energy, a 
more general rubric for cosmology in the early part of the 21st century might also 
be the "the distance revolution" . With new knowledge of the extra-galactic distance 
ladder we are, for the first time, beginning to accurately probe the cosmic expansion 
history beyond the local universe. While standard candles - most notably Type la 
supernovae (SNIa) - kicked off the revolution, it is clear that Statistical Standard 
Rulers, and the Baryon Acoustic Oscillations (BAO) in particular, will play an 
increasingly important role. 

In this review we cover the theoretical, observational and statistical aspects of the 
BAO as standard rulers and examine the impact BAO will have on our understand- 

t This review is an extended version of a chapter in the book Dark Energy Ed. Pilar Ruiz-Lapuente, 
Cambridge University Press (2010, ISBN-13: 9780521518888) 
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ing of dark energy, and the distance and expansion ladder. Fisher matrix forecasts 
for BAO surveys can be easily computed usingthe publically released, GUI-based, 
Fisher^Cast code, which is publically availabl4i|. 



1.1.1 A Brief History of Standard Rulers and the BA O 



Let us start by putting the BAO in context. The idea of a standard ruler is one 
familiar from everyday life. We judge the distance of an object of known length 
(such as a person) by its angular size. The further away it is, the smaller it appears. 
The same idea applies in cosmology, with one major complication: space can be 
curved. This is similar to trying to judge the distance of our known object through 
a smooth lens of unknown curvature. Now when it appears small, we are no longer 
sure it is because it is far away. It may be near and simply appear small because the 
lens is distorting the image. This degeneracy between the curvature of space and 
radial distance has not been the major practical complication in cosmology over the 
past century, however. That honour goes to a fact that has plagued us since the 
beginning of cosmology: we don't know how big extragalactic objects are in general, 
in the same way that we don't know how bright they intrinsically are. This problem 
was at the heart of the great debate between Shapley and Curtis over the nature 
of galaxies. Shapley argued that they were small and inside our own galaxy while 
Curtis maintained that they were extragalactic and hence much larger. 



To be useful for cosmology, we need a standard ruler, an object of a known size 
at a single redshift, z, or a population of objects at different redshifts whose size 
changes in a well-known way (or is actually constant) with redshift. Ideally the 
standard ruler falls into both classes, which, as we will argue below, is the case for 
the BAO, to good approximation. BAO are however a new addition to the family 
of putative standard rulers. A few that have been considered in the past include 
ultra-compact radio sources ({g^ . [s^ ) which indeed lead in 1996, prior to the SNIa 
results, to claims that the density of dark matter was low, 0^ < 0.3, with a non- 
zero cosmological constant of indeterminate sign (|67). Another radio standard ruler 
candidate is provided by double- lobed radio sources (j24l ). These Fanaroff- Riley Type 
lib radio galaxies were suggested as cosmological probes as early as 1994 (14 ih . and 
subsequent analyses have given results consistent with those from SNIa (j42l . |43| ). 
An alternative approach uses galaxy clusters. Allen et al. relate the X-ray flux to 
the cluster gas mass, and in turn, its size, providing another standard ruler under 
the assumption that the gas fraction is constant in time. This too leads to results 
consistent with those from SNIa (la, YII 



Beyond this we move into the realm of Statistical Standard Rulers (SSR), of which 
BAO are the archetype. SSR exploit the idea that the clustering of galaxies may 

I http://www.cosmology.org.za 
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have a preferred scale in it which, when observed at different redshifts, can be used 
to constrain the angular diameter distance. The idea of using a preferred clustering 
scale to learn about the expansion history of the cosmos has a fairly long history 
in cosmology, dating back at least to 1987 and perhaps earlier. In their conclusions 
Shanks et al. (jool ) forsee: 

There is one further important reason for searching for weak features in (r) , the quasar- 
quasar correlation function, at large separations. If a particular feature were found to appear 
in both the galaxy correlation function at low redshift and the QSO correlation function at 
high redshift, then a promising new cosmological test for qq might be possible. 



A series of later analyses further built up the idea of SSR in cosmology using vari- 
ously as motivation the turn-over in the power spectrum due to the transition from 
radiation to matter-domination, the mysterious 128/i~^Mpc feature detected in early 
pencil-beam surveys (|23 ) and the realisation that inflation could inject a preferred 
scale into the primordial power spectrum. These early studies often found tenta- 
tive evidence for a low-density universe and/or non-zero cosmological constant e.g. 
(j44l . I23I ) . Since the preferred scale could not be accurately predicted a priori these 
studies only provided the relative size of the SSR at different redshifts. Neverthe- 
less, it was realised that this could provide interesting constraints on the expansion 
history of the universe (89, [oT 



BAO entered the fray initially as a putative explanation for the apparent excess 
clustering around 100fe~^Mpc but were found to be too weak to be the origin for 
the apparent excess (|47|, l76|). The idea of using BAO themselves to learn about 
cosmological parameters seems to date first from Eisenstein et al. (1998) (j48l ) who 
wrote: 



Detection of acoustic oscillations in the matter power spectrum would be a triumph for cos- 
mology, as it would confirm the standard thermal history and the gravitational instability 
paradigm. Moreover, because the matter power spectrum displays these oscillations in a dif- 
ferent manner than does the CMB, we would gain new leverage on cosmological parameters. 



The first photometric proposal for using the BAO as standard rulers for learning 
about cosmology appears to date from 2001 (0). The real foundations, however, of 
the modern ideas on BAO, their detection and use, were laid by Eisenstein (2003) 
Blake and Glazebrook (2003) Jli) and a slew of later papers (jol, @, |2l|, [sO . 




lll|,[6i,0,[ 100 , 8^, Issl ) which developed hand-in-hand with the analysis of real 



data. Tantalising hints for the existence of the BAO were already visible in the 
Abell cluster catalogue (j77l ). but definitive detections had to wait for the increased 
survey volume and number density of galaxies achieved in the SDSS and 2dF redshift 
surveys, which immediatel y yie lded strong constraints on both curvature and dark 



energy at z < 0.5 (|52l . |32| . llOq . l83l . 1 54 ). Figures (jl.ip and (jl.2p show the original 
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evidence for the acoustic signature in the correlation function and power spectrum. 
Extracting the BAO scale from the matter power spectrum remains a thriving area 
of research in contemporary cosmology, as we discuss later in Section [1.51 on current 
and future BAO surveys. 




Comoving Separation (h ^ Mpc) 

Fig. 1.1. The Baryon Acoustic Peak (BAP) in the correlation function - the BAP is visible 
in the clustering of the SDSS LRG galaxy sample, and is sensitive to the matter density 
(shown are models with flmh^ ~ 0.12 (top), 0.13 (second) and 0.14 (third), all with 
fifc/i^ = 0.024). The bottom line without a BAP is the correlation function in the pure 
CDM model, whh ilt = 0. From Eisenstein et ai, 2005 IH). 



1.1.2 Cosmological Observables 

We now discuss the relevant cosmological observables that are derived from standard 
rulers in general, and the BAO in particular. The Baryon Acoustic Oscillations in 
the radial and tangential directions provide measurements of the Hubble parameter 
and angular diameter distance respectively. The Hubble parameter, H = a/a - 
where a is the scale factor of the universe - can be written in dimensionless form 
using the Friedmann equation as 

E{z) = ^ = yJVL^il + zf + J)de/(^) + + zY + n,,d(l + zY , (1.1) 
-no 

where f[z) is the dimensionless dark energy density and = — „^ 2 = 1 ~ (^m + 
^DE + ^rad) is the density parameter of curvature with $7^ = corresponding to 
a flat cosmos. Vtm,Vtrad a-re the matter and radiation densities with corresponding 
equations of state Wi = pi/ pi = 0, ^ for i = m, rad respectively. 
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Fig. 1.2. Baryon Acoustic Oscillations (BAO) in the SDSS power spectra - the BAP of the 
previous figure now becomes a series of oscillations in the matter power spectrum of the 
SDSS sample. The power spectrum is computed for both the main SDSS sample (bottom 
curve) and the LRG sample (top curve), illustrating how LRGs are significantly more 
biase d than average galaxies. The solid lines show the ACDM fits to the WMAP3 data 
(|104[ ). wh ile the dashed lines include nonlinear corrections. Figure from Tegmark ct ai, 

2006 (IToeh. 



If one treats the dark energy as a barotropic fluid with an equation of state with arbi- 
trary redshift dependence, w^z), the continuity equation can be directly integrated to 
give the evolution of the dimensionless dark energy density, f{z) = Pde/pde{z = 0), 
via 



exp 



l + z' 



(1.2) 



When we quote constraints on dark energy it will typically be in terms of the CPL 
parameterisation (2^, 71 ) 

2; 

w{z)=Wo + Wa— , (1.3) 



1 + Z 



which has 



(l + ^)3(l+«;o+«'a)gxp|-3u;„- 



+ z 



(1.4) 



Much of the quest of modern cosmology is to constrain the allowed range of w{z) (or 
f{z)) and hence use this to learn about physics beyond the standard model of parti- 
cle physics and General Relativity. Apart from direct measurements of the Hubble 
rate, one of the ways to constrain w{z) using cosmology is through distance mea- 
surements. Core to defining distances in the FLRW universe is the dimensionless, 
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radial, comoving distance: 

dz' 



X(-)-/ (1.5) 







One then builds the standard cosmological distances using xiz) to give the angular 
diameter distance, dA{z) via 

dA{z) = ———^-—== sin (y^-nkx{z)) (1-6) 

and the luminosity distance, dL{z), given via the distance duality as 

dL{z) = {l + z)^dA{z) (1.7) 



/\ 

dA{z), dL{z) 



Fig. 1.3. Curvature and xi^) define cosmological distances - In a flat Universe, the cosmo- 
logical distances are determined by x{z) oc dz' /E{z'). In a general FLRW model, however, 
spatial curvature bends the light rays away from straight lines and hence alters distances, 
meaning that one needs to know both Qk and xi^)- a result distance measurements 
always show a degeneracy between curvature (fJ/t) and dynamics {H{z)). 

The expression (Eq. (jl.6p ) for dA{z) holds for all values of the curvature, 0^, since 
for < the complex argument in Eq. (jl.6p converts the sin function to the sinh 
function. Hence the two key quantities that determine distances in cosmology are 
the dimensionless distance x{z) and 0^, shown schematically in Figure (jl.3p . The 
link (Eq. (|1.7p ) between dA{z) and diiz) holds in any metric theory of gravity as 
long as photon number is conser ved. This distance duality can be tested and used 
to look for exotic physics (11, 12, lid . 7^ . ^). 



Distances have a significant disadvantage over pure Hubble measurements: they 
require an integral over f{z) which is itself an integral over w{z). Hence, any inter- 
esting features in w{z) tend to be washed out in distance measurements. There is 
also another problem: if we look at Eq. (jl.6p for dA{z) we notice that if we make 
no assumptions about f{z), then even perfe ct distance measurements cannot break 



the degeneracy between f{z) and Q.^ (|ll4l ). This is not a fundamental problem 



if one assumes that w{z) has finite degrees of freedom, e.g. in Eq. (|1.3p . but one 
must remember that the degeneracy is being broken artificially by hand through 
ones choice of parameterisation and not by the data. As an example, consider Fig- 
ure (jl.4p . which shows two possible dark energy survey configurations; one with 1% 



Baryon Acoustic Oscillations 7 
measurements of the angular diameter distance at redshifts of z = 1,3, and another 
with measurements of both the Hubble parameter and angular diameter distance, 
but with double the errors on each observable. Even given the increase in the error 
on the observables, the dark energy constraints are significantly improved when in- 
cluding data from these complimentary probes when marginalising over curvature. 
In principle this degeneracy can be broken even with arbitrary f{z) by simultane- 




Fig. 1.4. Breaking the curvature-dark energy degeneracy - Error ellipses for the CPL pa- 
rameters wo,Wa, in two survey scenarios after marginalising over curvature: one consisting 
of two measurements at z = 1 and z = 3 of the angular diameter distance, with 1% errors on 
dA{z) (blue dashed curve), and another of measurements of H{z) and dA{z) (solid orange 
curve) where the errors have been doubled for both observables. In both surveys we assume 
a prior on curvature of 30, and Prior(J7m) = Prior(_ffo) ~ 1000. Even given a weak prior on 
curvature, combining measurements from multiple probes helps break the curvature-dark 
energy degeneracy. Figure produced using Fisher4Cast. 



ous measurements of both Hubble and distance; one can explicitly write 0^ in any 
FLRW model (with no recourse to the Einstein field equations) as (j3ll . l63l): 

lHoD(z)Y ^ ' 

where D is the dimensionless, transverse, comoving distance with D{z) = {c/Hq){1 + 
z)dA{z) . This relation gives the value of 0^ today, as a function of measurements 
at any redshift z. Hence this can be turned into a powerful test of the Copernican 
Principle (jsol ). Since flk is a single number, the right hand side will have the same 
value when measured at any redshift if one is in a FLRW background. If it is found 
to vary with redshift, then we do not live in a FLRW universe. 



The beauty of Baryon Acoustic Oscillations is that they provide both dA{z) and 
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H{z) using almost completely linear physics (unlike SNIa for example which involve 
highly complex, nonlinear, poorly understood stellar explosions). In addition, they 
offer the as yet unproven possibility of delivering constraints on growth though the 
change in the amplitude of the power spectrum. The time-dependence of the matter 
density perturbations, 5p/ p obeys the equation 

5 + 2H6 = AirGpmS. (1.9) 

The time-dependence of the growing mode of this equation is given by the growth 
function, G{z) which in a flat, ACDM model satisfies (j45l ): 

Gi.)^'^l-^.1±^. (1.10) 

while in a general universe with curvature and dark energy dynamics Eq. (jl.Op can 
be rewritten as (fl^): 

6" + ^ (l + ^ - w{x)nuE{x)^ ^ - l^U^)^ = 0, (1.11) 

in terms of the dimensionless scale factor x = a/ao = 1/(1 + z), where ao = 
c/HQ^y/TTk) is the radius of curvature of the universe. From Eq. (ILlOh we can 
see that the growth contains a mixture of Hubble and "distance" information - as a 
result measurements of growth are potentially powerful probes of dark energy. 



1.1.3 Statistical Standard Rulers 



To illustrate the idea underlying Statistical Standard Rulers (SSR), imagine that 
all galaxies were positioned at the intersections of a regular three-dimensional grid 
of known spacing L. Measuring angular diameter distances as a function of redshift 
would be trivial in this case (c.f. Eq. (jl.6p ) and we would also have the expansion 
rate as a function of redshift, measured at a discrete set of redshifts corresponding 
to the mid-points between galaxies. Now imagine that we start to randomly insert 
galaxies into this regular grid. As the number of randomly distributed galaxies 
increases the regular grid pattern will rapidly become hard to see by eye. However, 
the underlying grid pattern would still be detectable statistically, for example in the 
Fourier transform. 

However, a regular grid distributed throughout space would provide an absolute 
reference frame and would break the continuous homogeneity of space down to a 
discrete subgroup (formed by those translations which are multiples of the grid spac- 
ing L). To get to the core of SSR consider the following prescription for building up 
a galaxy distribution. Throw down a galaxy at random. Now with some fixed prob- 
ability, p, put another galaxy at a distance L (in any direction) from it. Using the 
new galaxy as starting point repeat this process until you have the desired number 
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of galaxies. Now there is no regular grid of galaxies but L is still a preferred length 
in the distribution of the galaxies, and forms an SSR. To reconstruct this preferred 
scale is a statistical problem. This is illustrated schematically in Figure (jl.Sp . which 
shows many rings of the same characteristic radius L superimposed on one another. 
This superposition of rings on the plane visually 'hides' the characteristic scale, as 
number of rings increases, and the sampling of each individual ring is reduced. 
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Fig. 1.5. Rings of power superposed. Schematic galaxy distribution formed by placing the 
galaxies on rings of the same characteristic radius L. The preferred radial scale is clearly 
visible in the left hand panel with many galaxies per ring. The right hand panel shows a 
more realistic scenario - with many rings and relatively few galaxies per ring, implying that 
the preferred scale can only be recovered statistically. 



BAO, as we discuss below, provide an elegant SSR hidden between the rest of 
the galaxy clustering, but they are not the only possible SSR's. Any preferred 
scale in the clustering provides either an opportunity to apply a relative (absolute) 
Alcock-Paczynski test (0) in the case the we don't (do) know a priori the size of 
the preferred scale, L. Other preferred scales include the Hubble scale at matter- 
radiation equality (which controls the scale of the turnover in the matter power 
spectrum) and the Silk damping scale. There may well be other preferred scales 
imprinted into the primordial clustering of matter. These can be naturally achieved 
if one inserts a short period of fast rolling into the otherwise slow-roll of inflation. 
The sudden change of inflaton velocity creates a bump in the matter power spectrum 
that can serve the same purpose as the BAO. The required fast-roll can be achieved 
in multi-field models of inflation. 

The SSR provided by the BAO has an additional advantage: it is primarily a linear 
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Fig. 1.6. The radial length of an object is given by cdz/H{z) where dz is the difference in 
redshift between the front and back of the object while the transverse size of the object is 
dA{z)9 and 9 is its angular size. If one knows that the object is spherical (but docs not 
know the actual diameter) then one has the Alcock-Paczynski test which gives the product 
dA{z)H{z) from measuring dz/9. If, as in the case of BAO, one can theoretically determine 
the diameter, one has the bonus of finding dA{z) and H{z) separately. 

physics phenomenon, which means we can ignore nonlinear effect^ to good approx- 
imation (we will discuss them later however) . This also means we can turn the BAO 
into a calibrated or absolute Alcock-Paczynski test since the characteristic scale of 
the BAO is set by the sound horizon at decoupling. As a result the angular diameter 
distance and Hubble rate can be obtained separately. The characteristic scale, s^^{z), 
along the line-of-sight provides a measurement of the Hubble parameter through 

H{z) = ^ (1.12) 
while the tangential mode provides a measurement of the angular diameter distance. 

This is illustrated in the schematic Figure. (jl.6p . The horizontal axis is Eq. (jl.l2p 
and the vertical axis is cAz/H(z), Eq. ()1.13p . 

While the AP test on its own constrains the product d^iz) x H{z), it is just one 
function, and so combining the measurements of cIa and H through the BAO provide 
tighter constraints on cosmological parameters. This is illustrated in Figure (jl.7p . 
which shows the error ellipse in the dark energy parameters from a hypothetical 



t These nonlinear effects include redshift space distortions and nonlinear gravitational clustering, which 
will in general change the spherical nature of the oscillation scale. 
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Fig. 1.7. Comparing the Alcock-Paczynski (AP) test and BAO constraints - on the dark 
energy parameters wq , Wa in the CPL paramctcrisation. The sohd blue and dashed dark 
brown lines elUpses illustrate the constraints from a pair of 1% measurements at z = 1, 3 of 
the angular diameter distance and Hubble parameter respectively. The filled beige ellipse 
shows the constraints from a BAO-likc survey, where both the Hubble rate and angular 
diameter distance are measured simultaneously at z = 1,3 (the errors on each have been 
increased by V2). In contrast, the dotted line illustrates constraints from the AP test: where 
1.5% measurements are made on the product H x dA &t the same redshifts. The AP test 
provides comparable constraints to using only a single measurement while the BAO provide 
the tightest constraints; they measure and H simultaneously. The assumed fiducial 
model is {Hq, Qm, Qk,wo, Wa) = (70, 0.3, 0, —1,0), and the priors on the model are given as 
Prior(Ho, f^m, nk,wo,Wa) = (10", 10", 10^, 0, 0). Figure produced using Fisher^Cast {H). 



galaxy redshift survey with constraints from the angular diameter distance and 
Hubble parameter and the product H x (Ia from the AP test. Constraints on the 
dark energy parameters from the AP test are similar to those from a single observ- 
able such as dyi, while the constraints are significantly improved when combining 
measurements of both H and cIa- 



One method of extracting a statistical scale from the clustering of galaxies is via 
the two-point correlation function, ^(r), which quantifies the excess clustering on 
a given scale relative to a uniform distribution with the same mean densi ty. T he 
correlation function of galaxies is approximately described by a power law ( lOSi ). 



(1.14) 



12 

with ro '-^ 5/i~^Mpc~^. 
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A characteristic scale in the clustering of galaxies will appear as a peak or dip 
in the correlation function, depending on whether there is an excess or deficiency 
of clustering at that scale. Any characteristic features will also be present in the 
power spectrum, since the correlation function and power spectrum (we consider 
for simplicity the simple 1-dimensional spherically averaged power spectrum) form 
a Fourier pair: 

/oo 
^{r) exp{—ikr)r'^dr . (1-15) 
-oo 

We will now see how features in the two functions are related. A 6 function at 
a characteristic scale, say r^,, in ^(r) will result in power spectrum oscillations, 
P{k) oc e"*'^^*, as can be seen in Figure (jl.Sp . These are the Baryon Acoustic 
Oscillations. 
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Fig. 1.8. Schematic illustration of the Fourier pairs ^(r), P{k). A sharp peak in the corre- 
lation function (left panel) corresponds to a series of oscillations in P{k) (right panel). The 
Baryon Acoustic Peak in the correlation function will induce characteristic Baryon Acoustic 
Oscillations in the power spectrum. 



These characteristic oscillations are powerful probes of dark energy. If we used only 
the radial Fourier modes we would obtain a measurement of H{z) while the purely 
transverse modes yield a pure measurement of dA{z). If we limited ourselves to just 
the purely radial and transverse modes we would be throwing away a large number 
of modes (and hence information) since most Fourier modes are not purely radial 
or transverse of course, but rather contain components of both. As a result, the 
measurements of H{z) and (i^(z) from actual BAO surveys are actually partially 
anti-correlated and hence not independent. This anti-correlation actually leads to 
stronger cosmological constraints than an uncorrelated analysis would suggest, as 
described in (jol). 

A nice example of their sensitivity and the issues involved in using BAO for pa- 
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Fig. 1.9. Changing cosmology moves the Baryon Acoustic Oscillations - the "wiggles-only" 
power spectra of the SDSS LRG survey. In each panel, a fiducial model with a different value 
of Vim has been used to convert the data from redshift to comoving distance, from ilm = 0.1 
(top) to rim = 0.4 (bottom). The solid lines in each panel show the CDM prediction for 
the BAO assuming the particular value of ri„i, while the dashed lines show the same model 
without the low-redshift small-scale damping term. In each case the baryon fraction is held 
fixed at 17%, and h — 0.73. As the matter density is changed, both the theoretical scale 
of the BAO changes, and the data moves around through the conversion from redshift to 
comoving distance. Figure from Percival et al., 2006 (j83f ). 



rameter estimation is evident in Fig. p.9|) from Percival et al. (|83l ) which shows 
both the SDSS LRG data and theoretical fits as a function of The panels each 
show the residual BAO oscillations relative to a non-oscillatory reference spectrum 
(a technique advanced in Blake et al. |l^ )). thus isolating the BAO information 
alone. The panels have a fixed baryon fraction of 17% relative to dark matter. An 
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important point to take away from this plot is that both the theory and the data 
vary with Q^n- 

This situation is different from more standard model fitting exercises where the data 
is static as the model parameters are varied. Here, to compute the data points, a 
redshift-distance relation (which depends on of course) must be assumed to move 
the points into physical, as opposed to redshift space, after which the power spectrum 
is constructed. In (0), 31 data sets and theoretical predictions were computed for 
various flat ACDM models with different values of i^m- The impact of varying 
is seen mostly clearly on the theoretical curves where the amplitude and frequency 
of oscillation is affected, both by the change to the underlying BAO scale (since 
the pb I ratio changes and hence the speed of sound prior to decoupling) and the 
redshift-distance relation which will move galaxies in the radial direction, bringing 
them either closer or moving them further away in distance. This dependence of the 
data on the assumed cosmology is typically ignored in parameter estimation studies 
using the BAO results, an a ppro ximation that is fairly good if one does not range 
far from the fiducial model ( lod . 82). The dilation resulting from using a different 



fiducial model would scale the A;-axis of the dimensionless power spectrum by the 
ratio (P. flOel): 

^ ^ dv{z) 

^fiducial V • ) 

1 /3 

where dv{z) = [d^iz]'^ cz / H (z)^ combines the radial and transverse dilation. 

Tegmark et al. calibrate this for the parameter range allowed by the combined 
SDSS-LRG and WMAP3 data, and find that corrections to the /c-scale are typically 
less than or equal to 3%. As an example, using a fiducial model with Vt^ = 0.25 
results in a bias of the measured value of $7™ of ~ 2%. 



1.1.4 Physics of the BAO 

Before recombination and decoupling the universe consisted of a hot plasma of pho- 
tons and baryons which were tightly coupled via Thomson scattering. The com- 
peting forces of radiation pressure and gravity set up oscillations in the photon 
fluid. If we consider a single, spherical density perturbation in the tightly coupled 
baryon-photon plasma it will propagate outwards as an acoustic wave with a speed 
Cs = c/-\/3(l -|- R), where R = 3pb/4:p^ oc Qb/i^ + z) (jsol ). At recombination the cos- 
mos becomes neutral and the pressure on the baryons is removed. The baryon wave 
stalls while the photons freely propagate away forming what we now observe as the 
Cosmic Microwave Background (CMB). The characteristic radius of the spherical 
shell formed when the baryon wave stalled is imprinted on the distribution of the 
baryons as a density excess. The baryons and dark matter interact though gravity. 
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Fig. 1.10. Snapshots of an evolving spherical density perturbation the radial 
mass profile as a function of comoving radius for an initially point-like overdensity located 
at the origin. The perturbations in the dark matter (black curve), baryons (blue curve), 
photons (red) and neutrinos (green) evolve from early times (z = 6824, top left) to long 
after decoupling {z = 10, bottom right). Initially the density perturbation propagates 
through the photons and baryons as a single pulse (top left-hand panel). The drag of the 
photons and baryons on the dark matter is visible in the top right panel; the dark matter 
only interacts gravitationally and therefore its perturbation lags behind that of the tightly 
coupled plasma. During recombination, however, the photons start to "leak" away from the 
baryons (middle left panel); and once recombination is complete (z = 470, middle right) 
the photons freely steam away leaving only a density perturbation in the baryons around 
ISOMpc, and a dark matter perturbation near the origin. In the bottom two panels we see 
the how the gravitational interaction between the dark matter and the baryons affects the 
peak: dark matter pulls the baryons to the peak in the density near zero radius, while the 
baryons continue to drag the dark matter overdensity towards the 150Mpc peak (bottom 
left), finally yielding a peak in the radial mass profile of the dark matter at the scale set 
by the distance the baryon-photon acoustic wave could have travelled in the time before 
recoupling. Figure taken from Eisenstein et al. (fSOh . 
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and so the dark matter also preferentially clumps on this scale. There is thus a in- 
creased probability that a galaxy will form somewhere in the higher density remains 
of the stalled baryon wave than either side of the shell. This scenario is illustrated 
in the iconic Figure (jl.lOp from Eisenstein et al. (fsol ). 



Had a galaxy formed at the centre of our initial density perturbation there would be 
a bump in the two-point correlation function at the radius s of our spherical shell, 
reflecting the higher probability of finding two galaxies separated by the distance s. 
The scale s is usually close to the sound horizorQ, the comoving distance a sound 
wave could have travelled in the photon-baryon fluid by the time of decoupling, and 
depends on the baryon and matter densities via 0): 



c^dz 



1 



2c 



H{z) y^v^ysi;^ 



In 



eq 



V 1 + Rrec + a/ Rrec + Rt 



eq 



(1.17) 



where R = 3ph/4:p^ oc VLi,h? / [l+z), z^q = i^m/^rad is the redshift of matter-radiation 
equality and "rec" refers to recombination. The CMB strongly constrains the matter 
and baryon densities at decoupling and hence the sound horizon, 146.8 it l.SMpc 
(ZSi Hence this scale is itself an excellent standard ruler as long as one can 



measure 17;, to high precision, and in the case where one allows for exotic radiation 
components in the early Universe, the redshift of equality (jsil). Of course, the early 
universe was permeated by many such spherical acoustic waves and hence the final 
density distribution is a linear superposition of the small-amplitude sound waves 
in the usual Green's function sense. Indeed this real space approach to cosmic 
perturbations based on Green's functions can be made rigorous, see Bashinsky and 
Bertschinger (fiol ). 



1.2 Forecasting BAO Constraints and Pov^er Spectrum Errors 



To make projections for the ability of a BAO survey to constrain cosmological pa- 
rameters, one can eith er simulate datagiven survey parameters or perform a Fisher 
matrix analysis, ( 105 . 96, [tI, 97, @, [hH, Isl). The Fisher matrix, Fij, defines how 
survey accuracy on physical observables (for e.g. the power spectrum or correlation 
function) translates into constraints on parameters of interest, 6i, such as curvature 
or dark energy, with the marginalised error on the i-th parameter of interest being 
given via the ii-component of the inverse of the Fisher matrix: 



A9i > ^{F^i 



(1.18) 



t This is perhaps a sUght misnomer since s is not the maximum distance that could have been travelled 

by any fluid. A simple scalar field has speed of sound equal to c and hence travels faster than the 

baryon-photon sound wave we are discussing here. 
X Confusion sometimes arises from quoting the BAO scale as both ~ ISOMpc and 105/i~^Mpc. Here 

h ~ 0.7 is the Hubble constant today in units of lOOkm/s/Mpc and provides the required conversion 

factor. 
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while in the case where one does not marginahse over any other parameters the best 
error is simply given by A6i > {Fii)~^/'^ which is always smaller. In the case where 
the likelihood is exactly Gaussian in the parameters then these error estimates are 
exact and not simply lower-bounds. 

In the case of the Baryon Acoustic Oscillations the Fisher^ matrix can be given 
directly in terms of derivatives of the power spectrum P[ k ) with respect to the 
cosmological parameters of interest, Oi (9^) as a function of the vector k: 

/•^™ ainP(T) 91nP(T)_,^ ,->, d^~k 
= L Vefr(fc) 



' f'^-' d[nP(k^d]nP(k^ 2rfrfMf, 

Assuming azimuthal symmetry for P along the line of sight means that P{k) can 
be described by ^ = k ■ f/k, with f is the unit vector along the line of sight and 
k = \k\. Eq. (jl.lOp is valid between wavenumbers between /cmin and a maximum 
wavenumber to exclude the non-linear regime. A conservative value for /cmax is 



0.1 /iMpc"^ at z = (15^). 



A crucial ingredient for the Fisher matrix above is the effective volume: 



nP{k, /i) 
nP{k,fi) + 1 



Purvey 5 (1.20) 



where the final expression assumes that the number density n(T^) is independent 
of position. The effective volume transforms the theoretical, survey-independent 
derivatives into specific predictions for the survey one is considering. The effective 
volume encodes the accuracy with which the survey under consideration can measure 
the power spectrum at different wavenumbers, k. To g ain i nsight into this key 



element, let us model the error on the power spectrum as (jlOSl . l2ll ) 



(5P 1 / 1 , , , 

1 + ^ , (1-21) 



P ^Jm \ nP 

where m is the total number of independent Fourier modes contributing to the 
measurement of the oscillation scale, and P = P{k*), k* ~ 0.2/iMpc~^ is the value 
of the power spectrum amplitude at an average scale k* , characteristic of the acoustic 
oscillations. From Eq. (jl.21|) we note that the two competing sources of error in 
reconstructing the baryon acoustic oscillation scale are cosmic variance and shot 
noise, represented by the first and second terms respectively. 
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1.2.1 Shot noise 



Our goal is to reconstruct the underlying dark matter distribution from discrete 
tracers such as galaxies. As is illustrated schematically in Figure (jl.lip . this re- 
construction is very difficult if we have few objects, an effect known as Poisson 
shot noise. Shot noise error decreases as the density of targets in a given volume 
increases: the complex underlying pattern in Figure (|1.11|) becomes clearer as the 
number of points increases and the pattern is sufficiently sampled. Increasing the 
number of targets requires longer integration times on the same patch of sky to go 
deeper, leading to a reduction in the area surveyed in a fixed observing time. 

Recently, an exciting proposal has been made to drastically reduce the shot noise 
term (jOa ) . In the context of the halo model, the power spectrum of dark matter is 
made up by summing the contributions from dark matter halos of different sizes. 
Under the assumption that galaxies form within the halos according to a given halo 
occupation distribution, one can derive a similar expression for the power spectrum 
of galaxies. The central idea in ([ol) is that nonlinear evolution of structure in the 
dark matter ensures that the power spectrum of dark matter behaves as k'^ for ^ 0, 
as opposed to the usual constant (fc") Poisson term. Using an N-body simulation 
with a scheme where central halo galaxies are weighted by halo mass, Seljak et al. 
IqsI ) find that the shot noise term is suppressed by a factor 10 — 30. 



1.2.2 Cosmic variance 

Cosmic variance is the error arising when we can't see the big picture: we cannot 
estimate the clustering on scales larger than our survey size. This is illustrated 
schematically in Figure (jl.l2p . where points are distributed according to a pat- 
tern consisting of Fourier modes with a variety of wavelengths and directions. The 
survey size increases in a clockwise direction, starting from top left, allowing the 
large-wavelength modes to become visible. To reconstruct these large-scale patterns 
requires an increased survey volume. In cosmology, it is impossible to keep increas- 
ing the size of the survey indefinitely as we are limited to the observable universe 
and so the values of the power spectrum on the very largest scales are fundamentally 
limited by this cosmic variance. To minimise the cosmic variance requires sampling 
the largest possible volume, which at fixed observing time implies spending as little 
time integrating on each field as possible, i.e. the exact opposite of the optimal 
strategy to minimise shot noise. The total error on the power spectrum is therefore 
a combination of the effects of the finite size of the survey and the number density 
of objects used to sample the underlying distribution. 

An interesting recent development is the suggestion by McDonald and Seljak (jrsl ) 
that by using multiple tracers of the underlying dark matter distribution, with differ- 
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Fig. 1.11. The effect of shot noise - as the number of galaxies in a survey increases, one is 
able to reconstruct the underlying pattern in the distribution of those points more reliably. 
This is illustrated in the progression from the top left hand panel (100 points) to the 
bottom right hand panel (100 000 points) which arc all drawn from the same probability 
distribution. As the number of points increases, the sub-structure of the pattern becomes 
visible. 



ent biases, one can constrain the redshift distortion parameter /3 = 6~^dln G/dln a, 
to a level of precision controlled only by the shot noise, with the exciting possibility 
that cosmic variance may not be the limiting factor for certain surveys. 



1.2.3 Redshift Errors 



To compute P{k), or ^(r) we need the redshifts of the objects in our survey. This can 
be achieved by taking a spectrum, which typically gives a highly accurate redshift, 
6z < 0.1%, or by using the colours of the object alone to get a photometric redshift, 
which typically gives projections of 5z ~ 3 — 5%, for 5 optical bands. The trade- 




Fig. 1.12. Illustrating cosmic variance with successive zooms - as the window expands 
(clockwise from top left), so large-scale modes in the distribution of points become visible, 
due to the reduction in cosmic variance. The number of points is the same in each panel 
while the volume increases by a factor of 4 x 10^ over the six panels. 



off when going from a spectroscopic to a purely photometric survey is to sacrifice 
redshift accuracy for greater depth, area and volume. How does degrading redshift 
affect the accuracy of the survey? Well, consider a fractional error az on the redshift 
z. From our discussion of the radial comoving distance and the Alcock-Paczynski 
effect, we know that an error 6z will result in an uncertainty in the radial position of 
6L = c6z/H(z). If we consider z ~ 1 we see that even if 6z/{l + z) = 0.01, the radial 
uncertainty is significant at around 1% of the Hubble scale or about 40/i~^Mpc. This 
is a large error when we are trying to measure the BAO on a scale of ~ 105/i~^Mpc. 
It is significantly larger than the systematic errors due to nonlinear effects and worse 
we cannot calibrate for it, as will probably be possible for the nonlinear effects. In 
addition one must include the fact that photometric redshift errors are typically 
non-Gaussian. They usually exhibit catastrophic wings at certain redshifts which 
reflect confusion with objects with similar colours at completely different redshifts, 
which means the likelihood of getting a "Sa" error is much larger than the Gaussian 
approximation would suggest. 

However, these redshift errors only affect our knowledge of the radial position of the 
galaxy. In the angular direction, astrometry errors (69) induce transverse distance 
errors of order 6L ~ 66dA{z). If we demand that this be less than l/i~^Mpc at 
all redshifts of potential interest {z < 4) we need astrometry accuracy better than 
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about 0.5' which is easily satisfied in all modern galaxy surveys. Hence even if our 
knowledge of the radial distribution of galaxies is significantly degraded in photo- 
metric surveys, the angular distribution is preserved. As a result measurements of 
the Hubble rate (derived from the radial BAO scale) in photo-z surveys are more 
strongly affected than measurements of cLa (derived from the transverse BAO scale). 

We can model the radial degradation as a suppression of the radial power spectrum 
by a factor exp[— (A;||crf.)^], where ctj. converts the redshift uncertainty into a physical 
distance as per our earlier discussion Gr = ^"^(l + z)/H{z). In this sense the effect 
of photo-z errors are similar to the nonlinear effects discussed later, which wipe out 
information about the higher order oscillations. This smearing is clearly visible in 
Fig. (jl.l3p which shows the progressive smoothing in the radial direction that occurs 
as the photometric redshift error increases. This smoothing smears out the baryon 
acoustic peak in the correlation function, as shown in Figure (jl.l4p . Simpson et 
al. suggest that in order to recover the dark energy equation of state w to 1%, the 



dispersion Uz must be known to within 10 (jlOll ) 



A simulation-calibrated fitting formula for the accuracy with which photo-z surveys 



can measure the power spectrum is given by Blake et al. (j2ll ). which is included in 
the FisherJ^Cast code (jl3l ). The key scaling is given by: 




(1.22) 



where V is the survey volume and Gr is, as before, the radial positional uncertainty 
due to the photo-z error. We see that sacrificing redshift accuracy has a similar effect 
to decreasing the survey volume. In particular, Blake and Bridle (2005) 0) estimate 
that with 5z/{l + z) = 0.03 a photometric galaxy survey requires approximately an 
order of magnitude greater area (and hence volume) to match the BAO accuracy 
achieved by a spectroscopic survey with the same magnitude limit (in the magnitude 
range r = 21—23). Recent claims suggest that there may be useful information in the 
small-scale power spectrum which reduces this to a factor five (j25l ) which would have 
important implications for the attractiveness of future photometric BAO surveys. 

In the limit of large photo-z errors, it is standard to bin the galaxies up into common 
redshift bins and project the galaxies in each bin onto a sphere at the central redshift 
of the bin. As a result one considers instead the angular power spectrum in each bin. 
This essentially throws away all of the radial information and is therefore maximally 



conservative ([18j, |7! 
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Fig. 1.13. The effect of redshift errors on the BAO. A 1/i^^Mpc thick slice of a simulated 
dark matter distribution is shown in real space (bottom panel) and in redshift space 
(second from bottom). The top two panels are the same distribution in redshift space, 
but with photometric redshift errors added to each galaxy. The second panel illustrates the 
effect of Gaussian errors of 0.3% (as expected for the Physics of the Accelerating Universe 
(PAU) survey), while the top panel has a 3% error. The redshift errors make it difficult to 
reconstruct the Baryon oscillation scale, which is shown in the bottom left hand corner of 
all panels (as a circle of radius 100/i~^Mpc). From Bemtez et al, 2008 (flSl ). 

1.3 Nonlinear theory 

Our analysis so far has been predicated on the belief that the BAO are reliable 
standard rulers, which is based on the fact that the BAO scale is ~ 105 /i^^Mpc 
which, within the context of FLRW models, is squarely in the linear regime since 
the quasi-linear regime only extends to about 30/i~^Mpc, even at z = (and is 
significantly smaller at higher redshift). However, every candidate standard ruler or 
candle has a limit beyond which it cannot be trusted, for either theoretical or obser- 
vational reasons. In the case of BAO the factors that contribute to the breakdown of 
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Fig. 1.14. The BAO signature smeared by photo-z errors - the reconstructed correlation 
function 'bump' at z = 0.5. The circles are the correlation function from a million sim- 
ulated LRG halos with mass M > 3.7 x 10^'^ Mq, from a simulated survey volume of 
V = 27/i~^Gpc'^. The reconstructed linear correlation function with 6 = 3 is shown by 
the dashed line, while the solid line shows the nonlinear prediction from renormalised 
perturbation theory (fSSl). The triangle, square and cross symbols show the measured 
correlation function after a Gaussian error of crz/(l + z) ~ 0.003,0.007 and 0.03 in the 
line-of-site direction is introduced. The corresponding solid lines are the analytical predic- 
tions for the damping of the Fourier space power spectrum from photometric errors using 
nonlinear corrections. Figure taken from Benftez et al, 2008 (flsh. 



confidence are nonlinear clustering and scale-dependent bias. But as we will discuss 
below, there are reasons to be optimistic even about these potential problems. 

A key advantage of the BAO as a cosmological probe is that nonlinearities such as 
those induced from nonlinear gravitational clustering induce predictable shifts in the 
oscillation scale and hence can be modelled both analytically and through numerical 
simulations. The effect of the nonlinearities can then be calibrated for, something 
which is not possible for many other standard rulers and candles. Here we briefly 
outline the effects of nonlinearity and techniques to correct for such nonlinearities. 



Different prescriptions exist for the method of using the Baryon Acoustic Oscilla- 
tions in the power spectrum as cosmological tools. The full Fourier space galaxy 
correlation method uses the entire power spectrum (including the shape) (0), but is 
sensitive to nonlinearities such as scale-dependent bias and nonlinear redshift space 
distortions. The effect of these systematics is reduced if one removes the overall 



shape of the power spectrum by dividing by some reference cosmology ([211, l98l ) 



however you also lose any information contained in the overall shape and amplitude 
and so constraints on cosmological parameters will be weaker. 
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1.3.1 Nonlinear bias 



Measurements of galaxy clustering from redshift surveys yield the galaxy power spec- 
trum, which is traditionally related to the power spectrum of dark matter Pdm(^)) 
(which we are interested in) through the bias z), 

P^.i{k, z) = b\k, z)PBM{k, z) , (1.23) 

which in principle can be both redshift and scale dependent (53, 33, 113l . [g^ . lod . 



9l|,l37|). 



Even a moderate scale dependent bias will shift the peaks of the BAO and cause a 
systematic error in the standard ruler. In the extreme case one could even imprint 
oscillations not present in the underlying dark matter distribution. Fortunately 
there is a way out of this degeneracy. Clustering in redshift space is anisotropic 
due to redshift distortions. The radial component of the galaxy peculiar velocity 
contaminates the cosmological redshift in a characteristic, scale-dependent manner, 
which means that the power spectrum in redshift space is not isotropic, P(k) ^ P{k). 
On large scales galaxies falling into overdensities (such as clusters) are 'squashed' 
along the line sight (the Kaiser effect), while on scales smaller than clusters the 
velocity dispersion of the galaxies within the cluster leads to the 'finger of god' 
effect - clusters appear elongated along the line of sight 



We can expand the anisotropic power spectrum as 

P(k,z)= Plik,z)^l{f^) (1-24) 

Z=0,2,4,... 

where £;(//) are the Legendre polynomials, fj, = cos{9),k = |k| and the monopole 
Po{k,z) is the spherically averaged power spectrum we have been discussing for 
most of this review. The odd moments vanish by symmetry. Studies have shown 
that the extra information in the higher order moments Pi allow the recovery of 
essentially all the standard ruler information, even marginalising over a reasonable 
redshift an d sc ale-dependent bias (e.g. a four-parameter model), with future ex- 



periments (|ll7l ). To understand why the different multipoles would break the bias 
degeneracy, remember that the amplitude of the redshift distortions is controlled 
by the parameter (3 = Jlm/fc, where 7 ~ 0.6. Imagine an observed monopole galaxy 
power spectrum. If 6 — > then amplitude of the dark matter power spectrum must 
increase to leave the galaxy clustering unchanged. The larger dark matter cluster- 
ing will lead to larger velocities and hence larger redshift distortions which will be 
visible in the dipole and quadrupole power spectra. Including information from the 
full p ower spectrum allows one to calibrate for such a scale-dependent galaxy bias 



(Ilia l84l^. 



Before moving on we note that the dark matter power spectrum is the product of 
the initial power spectrum of the Universe, the growth function G{z) (defined in 
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Eq. (|1.1U|) ) and the transfer function T{k): 

P^uik, z) = G\z)T\k)Pi{k) . (1.25) 

It is clear from Eq. (jl.25p that at the level of the monopole power spectrum the 
growth function is completely degenerate with the a general bias. Redshift distor- 
tions and non-Gaussian clustering (measured e.g. through the bispectrum) offer 
the opportunity of determining the growth independent of the bias using the same 
principle we have discussed above. 



1.3.2 Movement and broadening of the peak 

The main systematic error due to nonlinearity is the shift in the peak of the 
corr elati o n fu nction due to mode-mode coupling, as has been studied extensively 
Toi, [io3). There are a couple of effects at play here. First, if the broadband 



correlation function (i.e. the smooth part without the peak) changes with time, the 
acoustic peak will shift too, simply due to elementary calculus. Secondly, consider 
the simple physical model introduced in Section [1.1. 41 where we thought of the cor- 
relation function peak as arising from an acoustic wave that moves outwards before 
stalling at recombination with a galaxy at the origin from which the spherical shell 
expanded. This relies on the insight that the two-point correlation function, which 
is a joint probability, can be rewritten as a conditional probability: what is the 
probability of finding a galaxy at distance, s, given that there is a galaxy at the 
origin. If we then consider the nonlinear evolution of such a sharp density shell, at 
rest shortly after recombination, we would expect it to undergo some collapse over 
the history of the cosmos due to its own self-gravity and due to the gravity provided 
by the galaxy at the centre thereby shrinking the radius of the shell and hence the 
standard ruler length, by a small but systematic amount of 1 — 3%. 

As can be seen in Fig. (jl.lSp . nonlinearities not only shift the peak, but also smooth 
out and broaden the peak of the correlation function. We can again understand this 
qualitatively using our simple model. Imagine a galaxy forms on the BAP shell. 
At the nonlinear level this galaxy is subject to the combined interactions with all 
other over and under-densities. In any given density realisation this net force may 
pull the galaxy outwards or inwards. Since the correlation function is computed by 
averaging over all available galaxies, the average effect is to broaden the BAP. The 
only constant in all the shells is the galaxy at the centre which causes the small 
inwards shift of the BAP as described above. 

This broadening of the BAP equivalently can be thought of damping the oscillations 
in the power spectrurr[3 on small scales. Broadening the peak obviously makes 
reconstruction of the position of the peak - and hence the standard ruler length - 

t The correlation function and power spectrum arc introduced in Section ll.1.31 
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Fig. 1.15. Movement of the Baryon Acoustic Peak. The correlation function at z=0 illus- 
trates how the transfer of power to smaller scales due to nonlinearities, leads to a shift of 
the BAP in the correlation function, £,{r). The linear peak, indicated by the dashed line is 
both broadened and shifted towards smaller scales. The solid line shows the prediction for 
the shift from renormalised perturbation theory (RPT) (fSSl . [39t). The vertical lines denote 
the corresponding maxima of the linear and nonlinear correlation functions. From Crocce 
et ai, 2008 g^). 



less accurate, hence degrading dark energy constraints. Broadening the peak in the 
correlation function washes out the oscillations in -P(fc) at large wavenumbers or 
small scales. 



We can illustrate this analytically as follows. Let us model the correlation function 
as a Gaussian bump shifted so it is centered at a scale r*, or 

e(r)=expf-i!l^V (1.26) 



Hence the power spectrum is given by 



oo 



(^ — ^* ) ^ 

P{k) = I exp ( ^—^ — )exp{—ikr)dr 



oo 



= .^^^exp (— iA;r*) exp (— fc^iT^) (1-27) 

Figure (jl.lOp illustrates this toy-model correlation function consisting of a Gaussian 
shifted to some preferred scale r^,, and the corresponding power spectrum P{k). 
The oscillations are given for a range of widths of the Gaussian bump, 10 < ex < 35. 
Clearly as the Gaussian broadens, the oscillations in the power spectrum are washed 
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out, making their detection harder. Recalhng the illustration of rings of power in 
Figure (jl.Sp . we can examine the effect of successively broadening the rings from 
which the points are drawn. This this shown in Figure (jl.lTp which shows the 
smearing of the characteristic radius implying an increased error in the standard 
ruler measurement. 




Fig. 1.16. Smoothing out the baryon acoustic signal in the Fourier pair ^(r), P{k). Increasing 
the width of the correlation function bump corresponds to the damping of the acoustic 
oscillations in the power spectrum, particularly severely at large k. Both clfccts make 
reconstruction of the standard ruler length more noisy. 



1.3.3 Reconstruction 

While non-linear gravitational collapse broadens and shifts the peak of the corre- 
lation function, Eisenstein et al. (|4£|) point out that the map of galaxies used to 
extract the power spectrum in redshift space can also be used to map the veloc- 
ity field. Since the galaxies are essentially test particles in the standard ACDM 
paradigm, this velocity field can then be used to undo the effects of the nonlinear 
clustering or equivalently to reconstruct the position and sharpness of the linear 
acoustic oscillation peak by moving densities to where they would have been had 
linear theory held at all times. By considering a pair of galaxies separated by the 
characteristic BAO scale, Eisenstein, Seo and White (2007) (jsol ) show that the ma- 
jority of the corrupting signal comes from wavenumbers k ~ 0.02 — 0.2/iMpc~^. 
Larger wavelengths coherently move both galaxies while smaller scales are weak be- 
cause the power spectrum has little power there. The typical distances induced by 
nonlinear corrections are around 10/i~^Mpc. 

Various methods can be followed to reconstruct the velocity and density field which 
are summarised in (|4^). Eisenstein et al. move the measured densities back to their 
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Fig. 1.17. Hiding the characteristic scale. As the peak is broadened from top to bottom as 
shown schematically in Figure (|1.16p . the underlying rings of power are lost, and must be 
recovered statistically, c.f. Figure (|1.5p . The number of points are kept the same in each 
panel. 
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Fig. 1.18. Reconstructing the Baryon Acoustic Peak (BAP) - The nonlinear effects generate 
bulk flows that can be reconstructed from the galaxy distribution itself hence allowing the 
nonlinear movements of densities to be undone to good accuracy, thereby both sharpening 
the acoustic peak and moving it to the correct linear position. From Eisenstein et at, 2007 



linear locations using the following prescription. First they smooth the density field 
on about 10/i~^Mpc scales. Then they compute the Lagrangian displacement field, 
~q which is assumed to be irrotational (no vector perturbations) and which obeys 
the linear prediction A - ~q = —S. All particles are then shifted by — In redshift 
space the densities are then boosted by 1 + d{ln G) / d{ln a) , where G is the growth 
factor, to account for the linear redshift distortions. 

Figure (jl.lSp shows the reconstruction of the acoustic peak in the real-space corre- 
lation function using these techniques with improvement in the accuracy of the peak 
by a factor of 2 to 3, with similar results for the redshift-space correlation function. 

The impact of reconstruction on BAO survey optimisation is still to be explored 
but an immediate implication is that low-redshift surveys, where the nonlinearities 
and broadening of the peak is stronger, will benefit more from reconstruction than 
high-redshift surveys where the effect of nonlinearities have not had time to imprint 
on the relevant scales. This will, in turn, make low-redshift BAO surveys more 
interesting, allowing smaller telescopes (e.g. 2-4m class) to compete with the larger 
telescopes (8-lOm class) in the BAO stakes. We now move on to a discussion of 
targets for BAO surveys. 



61). 
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1.4 Target Selection 
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Fig. 1.19. Exposure time required to reach a given galaxy density on an 8m class telescope 
for four typical BAO targets: red galaxies at z ~ 1 (solid line) where the continuum is used 
to obtain redshifts; blue galaxies at z ~ 1 (dashed line) where redshifts arc obtained using 
galaxy emission lines; red galaxies at high redshift (z ~ 3 - dot-dashed line) and blue 
galaxies again at high redshift (dotted line). The plateau in the galaxy density occurs 
when the surface density reaches the spectroscopic fibre density (a single pointing of the 
telescope is assumed). Figure from Parkinson et ai, 2006 (|80l ) for the WFMOS survey. 



A key decision in undertaking any BAO survey is the choice of target, since the 
bias, b, of different potential targets differ considerably as a function of morphology, 
colour etc... relative to the underlying dark matter distribution, in addition to any 
redshift or scale dependence; Ptarget{k) = b'^PoMik)- This translates into different 
optimal target densities since if one requires nPtarget ~ 1 (the criterion translates 
to n ~ 1/(6^Pdm)): one needs fewer more highly biased tracers of clustering than 
weakly biased targets. For example. Luminous Red Galaxies (LRGs) are highly 
biased tracers, Iilrg ^ 

1.5 - 2 (Hoi) since they are typically found in clusters while 
blue spirals are typically field galaxies and hence are not strongly biased. All modern 
studies of BAO use dedicated targets, the choice of which typically trades off bias 
versus integration time. Integration times for various types of possible targets for a 
large 5000 fibre 10m class survey (such as WFMOS) are shown in Fig. p.l9p . 



1.4^.1 Luminous Red Galaxies (LRGs) 



Luminous Red Galaxies rose to prominence with the SDSS LRG survey (j52l ). They 
are typically "red and dead," passive elliptical galaxies with featureless spectra. A 
high S/N LRG spectrum is shown in Figure (jl.20p . The redshift is derived from the 
position of the 4000 A break which therefore requires long integration times even on 
a 10m class telescope, for redshifts z > 1. This is counteracted by the large bias of 
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RA=324,52667, DEC= 0,27758, UJD=52468, PlotB= 989, Fiber=596 




Wavelength [1] 

Fig. 1.20. Example of a Luminous Red Galaxy. The Sloan Digital Sky Survey DR6 image 
(top panel) and spectrum (bottom panel) of a typical LRG and showing a lack of 
emission lines, which means that the continuum must be measured to obtain a redshift. 



Prom http: / / skyserver.sdss.org/| 



LRGs which means that the required target density is significantly lower. The latter 
advantage, plus the ability to efficiently find LRGs in optical photometric surveys 
like the SDSS survey has lead to LRGs been chosen as the targets for the BOSS 
SDSS-II Survey 



1.4-2 Blue galaxies 

(|2ol). While redshifts for LRGs are obtained from the continuum spectrum, blue, 
star-forming galaxies have strong emission lines which provide good redshifts (see 
Fig. (|1.2ip ) from for example the Oil doublet at 3727 A which is within the optical 
band at redshifts z < 1.4 - expected integration times for this line using a 10m 
class telescope are around 15 minutes. Despite this, the low bias means that a 
much higher target density is required compared to LRGs. Selection of star- forming 
targets is achieved with a combination of optical and UV imaging and forms the 
basis for the WiggleZ survey which uses a combination of SDSS and GALEX (UV) 
imaging for selection. 



j www.sdss3.org 

X http: / / sdss3.org/ collaboration / description.pdf 
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4000 4500 5000 

Wavelength (Angstrom) 

Fig. 1.21. A star-forming galaxy spectrum showing the characteristic emission lines used for 
redshift determination taken as part of the WiggleZ survey l(20t ) . 



1.4-3 Lyman Break Galaxies 



The standard emission lines go out of the optical band at z ~ 1.4 leadmg to the 
redshift desert for optical surveys because of a dearth of emission lines at wavelengths 
< 3000A . This "drought" is broken by Ly-a at the wavelength of 1216A which 
comes into the optical passbands around a redshift of 2: ~ 2.3 and remains there until 
z ~ 6.4, making it an ideal target at high redshift. For galaxies at higher redshifts 
the Ly— a break moves into different bands and the galaxy will have negligible flux 
in (for example) the U band, but strong flux in the V band - hence the UV 'drop- 
out', multi-colour imaging of the galaxy is hence used to photometrically determine 
the redshift of the galaxy. Lyman Break Galaxies (LBG) take long integration times 
(see the high-z continuum curve in Fig. (|1.19p ) but there are large numbers of them. 



1.4 -4 Lyman Emitting Galaxies 



A small set of LBGs also have strong Ly-a emission lines. When they exist they 
provide ideal targets for redshifts due to the strong emission, however their number 
density is somewhat unknown. They are the target of preference for the HETDEX^ 
({giI ) instrument planned for the Hobby-Eberly Telescope in the redshift interval 
1.9 < z < 3.5, which is expected to detect 8 x 10^ of these Ly-a emitting galaxies. 



t http: / /hotdex.org 
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LSST will potentially detect millions of photometric Type la supernovae (SNIa). 
Zhan et al. have proposed that these supernovae could be used to measure BAO 
( II9I ) at z < 1. Hence the same data could provide both d/, and dA, providing 
constraints that are insensitive to Cosmic Microwave Background priors. The ad- 
vantage of SNIa as photometric BAO targets over galaxies is that the photometric 
redshift error is typically significantly smaller [gz ~ 0.02) due to the well-sampled, 
multi-epoch spectral templates that will be available from current and future low-z 
supernova surveys. 



1.4 -6 Lyman Alpha Forest 



Sampling the underlying dark matter distribution at a discrete set of N points 
makes it difficult to uncover subtle underlying patterns due to shot noise. Instead, 
a potentially superior method would be to take 1-dimensional slices through the 
density distribution. This is the idea behind using the Lyman-alpha forest to probe 
the BAO. McDonald and Eisenstein (2007) ({t^ ) discuss such a survey at redshift 
2.2 < z < 3.3, and project constraints on the radial and tangential oscillation scales 
of order 1.4%. Such a survey could be performed at the same time as a spectroscopic 
galaxy redshift survey, which is the plan for the BOSS SDSS-III survey. 



1.4 '7 21cm Neutral Hydrogen 

One step better than a 1-D slice is a full 3-D slice through the neutral hydrogen 
distribution. This is the possibility affo rded by neutral hydrogen surveys based on 
the 21cm HI emission line (fisl . [3, I1I6 . 27). The advantage of this probe is that 



neutral hydrogen should be ubiquitous at all redshifts, although the precise redshift 
dependence of the HI density is unknown and is further affected by uncertainties in 
cosmic reionisation. However, in principle, HI surveys will be able to probe deep 
into the dark ages before the formation of galaxies, providing access to the cosmic 
density field uncontaminated by nonlinearities. 



1.5 Current and Future BAO Surveys 



The key elements for a BAO survey are redshift accuracy, redshift coverage, area 
and volume (of course the latter three are not independent at fixed total survey 
time) . The ideal instrument therefore has large field of view (the area it can see at 
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any one time), large mirror size allowing short integration times and if it is taking 
spectra, the ability to harvest large numbers of spectra simultaneously. For future 
surveys being considered now, this means fields of view in excess of 1 deg^, mirrors 
greater than 4m in size and for spectroscopic surveys, the ability to take at least 
1000 spectra simultaneously, using a multi-fibre or other technology. 



1.5.1 Spectroscopic Surveys 

We now discuss in rough chronological order the current and future spectroscopic 
BAO survey^. First up are the final Sloan Digital Sky Survey (SDSS-II) LRG and 
main galaxy surveys at z ~ 0.35 and 2; ~ 0.1 respectively and which will cover 
about 10, 000 deg^ in the northern hemisphere. Next is the WiggleZ survey us- 
ing the 400 fibres on the AAT (|2ol ) and covering 1000 deg^ over the redshift range 
0.2 < z < 1.0, with a median redshift of 2; = 0.6, which will be completed in 2010 
and will measure H{z) and (iyi(z) to around 5%. Beyond that is the Baryon Oscil- 
lation Spectroscopic Survey (BOSS), part of the SDSS-IIl[il survey and scheduled 
to operate over the period 2009-2014. The relatively small diameter (2.5m) of the 
SDSS telescope combined with the large field of view means that BOSS will focus 
on a wide-and-shallow survey measuring approximately 1.5 million LRGs at z < 0.7 
and around 160,000 Ly-a forest lines at 2.2 < z < 3 and giving projected absolute 
distance measurements of 1.0% at z = 0.35, 1.1% at z = 0.6, and 1.5% at z = 2.5 [§|. 

Projects on a similar, 2010-2015, time-scale include FMOS and LAMOST. FMOS 
is an infra-red spectrograph for Subaru with 400 fibers which could undertake a 
moderate but interesting BAO survey in the redshift desert at z ~ 1 — 1.7 over 
~ 300 deg^ while LAMOST is the Chinese 4m telescope with a 4000-fiber spec- 
troscopic and 20 deg^ field of view which should enable a v ery e ffective BAO survey 
similar to BOSS both at 2; ~ 1 and high-z using quasar^ (|ll2l ). 

Another exciting BAO survey is the Hobby-Eberly Dark Energy eXperiment (HET- 
(j6ll ) which will target the highly biased Ly-a emitting galaxies over the 
range 1.8 < z < 3.7. Such a survey over 200 deg^ would probe about 5h~'^Gpc^ 
with approximately one million galaxies, allowing HETDEX to provide ~ 1% mea- 
surements of dA{z) and H{z) at three redshifts over the survey range. An attractive 
feature of HETDEX is that it does not need any pre-selection imaging; targets are 
acquired purely by chance using integral field spectrographs (jgl). 

t These would be called Stage II, III and IV surveys in the DETF report Q). 

X www.sdss3.org 

§ http:/ /sdssS.org/eoUaboration/description.pdf 

*i http://www.sstd.rl.ac.uk/fmos/ 

II http://www.lamost.org/en/ 

tf http: / /hetdex.org 
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The next major advance in the spectroscopic BAO domain would be enabled by the 
Wide Field Multi-Object Spectrograph (WFMOS) on a 10m class telescope, such as 



Subaru (|57l . Il4l). While WFMOS has been cancelled as a Gemini project, it may still 
take place in a different form. The default plan for WFMOS called for a large field 
of view (> 1 deg^) and a large number of fibers (at least 3000, although the optimal 
number is still being investigated (jsol )). Slated for a ~ 2015 start, the default 
WFMOS-like surveys would measure in excess of one million LRGs or blue galaxies 
at z = 0.8 — 1.3 over an area of 2000 - 6000 deg^ and of order one million Lyman 
Break Galaxies at z = 2.5 — 4 over a somewhat smaller area, providin g pe rcent level 



measurements of both dA{z) and H{z) at z = 1 and z = 3 (|46l . |9g, |80|, l8ll ). The 
high-z component of the survey would provide a key leverage against uncertainties 
in curva ture and w{z) at z > 1. It would also be a powerful probe of modified 



gravity (|ll7l ) and allow high spectral resolution archeology of the Milky Way to 



understand the origins of its stellar populations (|57j). 

Beyond the 2015 timescale there are a number of planned and proposed missions 
in various stages of preparation. Perhaps the simplest proposal is the BigBOSS 
ground-based experiment which, over a ten-year period would be able to compete 



with Dark Energy Task Force Stage IV experiments (j94l ). Building on BOSS it 
would use 4m telescopes at both northern and southern NOAO sites (initially 6 
years at KPNO potentially followed by four years at CTIO after the Dark Energy 
Survey) fitted with a new 4000-fiber R=5000 spectrograph covering 340 — 1130nm 
with a 7 deg^ field of view, yielding a survey sample of up to 50 million galaxies and 
a million quasars over 24000 deg^ and allowing superb measurements of the BAO 
and redshift space distortions in the range 0.2 < z < 3.5. At low-z, targets would be 
LRGs while for 1 < z < 2 BigBOSS would target bright Oil emission line galaxies 
with the QSOs taking over at z > 2. 



A more radical proposal is to use slitless spectroscopy (j55l . 1871 ) which is one possi- 
bility for the spectroscopy component of the proposed EUCLID survey, which is a 
combination of the earlier SPACE (0) and DUNE (jsol ) missions. The spectroscopic 
component posits an all-sky near-IR survey down to H=22 which would provide 
of order 150 million redshifts. EUCLID would aim for launch around 2018 if it is 
chosen as the winner of the ESA Cosmic Visions program. A further space BAO 
proposal for the DOE-NASA JDEM mission is ADEPT which would also gather 
around 100 million redshifts over the redshift range z < 2. Recently, the possibility 
of a JDEM-EUCLID merger has been raised due to the obvious complementarity of 
the science of the two programs and budgetary constraints, although the technical 
and organisational challenges of building such a complex joint US-Europe mission 
are likely to be significant. 



A very different direction is provided by radio BAO surveys. Despite the inherent 
weakness of the 21cm signal it is likely that radio telescopes will play an important 
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role in future cosmology. This is primarily driven by the fact that the sensitivity of 
radio telescopes for projects such as BAO scales as the square of their area, unlike 
optical telescopes whose sensitivity scale linearly with their diameter. This, together 
with technologies such as synthetic aperture arrays that allow very large fields of 
view offer the appealing possibility of surveying huge volumes at very high target 
densities. 



An exciting proposal in this direction is the 21cm Hubble Sphere Hydrogen Survey 
(HSHS) 3 which would measure the BAO in neutral hydrogen over the whole sky 
out to z = 1.5. This highly ambitious proposal would provide essentially cosmic 
variance- limited measurements of the power spectrum in bins of width Az ~ 0.1 
and exquisite accuracy on and H{z) in the same bins. The key to the HSHS 
concept is simultaneously combining huge collecting area with very a large field of 
view. This can be achieved, at what is hoped to be low cost, by using multiple 
fixed parabolic cylinders which provide drift scans of the entire sky everyday. In 
this sense, one of the Fourier transforms needed to form an image is undertaken in 
software ('along the cylinder') while the other is done in hardware ('in the parabolic 
direction'). HSHS is unusual for a galaxy survey because of its low angular resolution 
of around 1', adapted for statistical analysis of the BAO rather than producing a 
galaxy catalogue as its primary output. In this sense HSHS resembles a CMB 
experiment for neutral hydrogen. 

A more ambitious proposal is that of the full Square Kilometer Array (SKA)[3 which 
may be a fully software telescope at 21cm frequencies, with both Fourier transforms 
being done in software and using completely flat reflectors. The great advantage of 
such purely synthetic apertures would be that detectors would essentially see all of 
the visible sky all of the time, providing the ultimate field of view (0). This idea ap- 
pears to have been rediscovered in the form of the Fast Fourier Transform Telescope 



(|l07l ). The SKA would provide essentially cosmic variance limited BAO measure- 
ments out to z = 1.4 and beyond with of order 10^ redshifts, but with sub-arcsecond 
angular resolution, allowing in addition excellent weak lensing measurements (fl6l ). 
While SKA will be an exceptional BAO machine, pathfinders leadirig up to the full 
SKA will also provide the first detections of the BAO in the radio yj) . 



Beyond SKA one can imagine using radio surveys to probe the BAO at very high 
redshifts, z > 10, where many more modes are in the linear regime. Since there are 
essentially no galaxies above this redshift, neutral hydrogen will likely be the only 
way to test dark energy in the dark ages at (jisl). 



t http://hlsurvcy.phys.cmu.edu/ 
I http://www.skatcloscopo.org/ 
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Spectra are slow and expensive to obtain and it is tempting to try to study the 
BAO with only multi-band imaging. A large number of photometrically harvested 
galaxies might provide a useful probe of the Baryon Acoustic Oscillations, provided 
the photometric redshift error is small enough, as discussed in Section 11.2.31 



The current state of the art of pho tometric redshift surveys is provided by the MegaZ 
and related catalogues ( 34. 1 181.1 73) based on the SDSS photometry. These catalogues 
target LRGs and typically achieve 5z ~ 0.03(1 + z) redshift accuracy with approx- 
imately 1% contamination from M-star interlopers after suitable cuts. Although 
they include more than 1 million LRGs out to z ~ 0.7 and cover 10,000 deg^ they 
do not detect the BAO with any significance due to projection effects arising from 
the photometric redshift errors. 



Beyond SDSS there are a number of exciting photometric surveys. SkyMapper will 
essentially provide Sloan in the southern hemisphere (joil ) while the Dark Energy 
Survey (DES) @) will use 30% of the 4m CTIO telescope time to cover around 
5000 deg^ and detect of order 300 million galaxies in the five Sloan photometric 
bands, u,g,r,i,z over the redshift range 0.2 < z < 1.3 while the PSl phase of 
the Pan-STARRS project could cover Stt steradians of the sky and detect of order 
100 million LRGs, again in five, slightly redder, passbands (|25l). Both surveys 
should provide compelling BAO detections in addition to the wealth of other science 
including lensing and a rich SNIa dataset. 



A further interesting hybrid is the Physics of the Accelerating Universe (PAU) pro- 
posal which one might call an ultra-photometric or quasi-spectroscopic survey. PAU 
plans to bridge the gap between standard photometric and spectroscopic surveys 
through the use of an order of magnitude more filters than SDSS, DES or Pan- 
STARRS. Using 40 narrow-band and two broad-band filters covering the optical 
range, the aim is to identify the 4000A break with enough spectral resolution to 
determine redshifts to an accuracy of 6z ~ 0.003(1 + z) which is hoped will provide 
sufficient accuracy to reconstruct the BAO scale in the radial direction and hence 
obtain H{z) information as well as dAiz) (|l5l ). Again LRGs are the targets of choice 
due to their simple spectra and with a survey area of order 8000 deg^ the desire is 
to measure such ultra-photometric redshifts for over 10'' LRGs at z < 1, although 
doubts have been raised as to whether this approach is competitive with spectro- 
scopic BAO surveys (jssl ). The correlation function from simulated galaxy halos is 
given in Figure (|1.14p . from (fisl). The smearing of the acoustic feature is clearly 
visible as the photometric redshift error increases. 



Beyond these surveys the Large Synoptic Survey Telescope (LSST) will likely provide 



t https:/ /www. darkenergysurvey.org/ 
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the definitive photometric survey for the next two decades. Covering 20, 000 deg^ 
of the sky visible from Chile, LSST would detect every galaxy visible in the optical 
down to a co-added limiting magnitude of r = 27.5, or about 10 billion galaxies. 
With science operations slated to begin in 2015 or soon thereafter, LSST will yield 
exquisite detections of the angular BAO as a function of redshift, albeit wit hout the 
radial information provided by spectroscopic or ultra-photometric surveys (Hoi). 



1.6 Conclusions 



In the era of precision cosmology, standard rulers of ever-increasing accuracy will 
provide powerful constraints on dark energy and other cosmic parameters. The 
Baryon Acoustic Oscillations are rooted primarily in linear physics with nonlineari- 
ties that can be well-modelled and corrected for. As a result the characteristic scale 
of these 'frozen relics' imprinted into the cosmic plasma before decoupling will likely 
remain as the most reliable of the Statistical Standard Rulers in the coming decade. 
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